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Abstract 

We discuss flavor-mixing probabilities and flavor ratios of high energy astrophysical 
neutrinos. In the first part of this paper, we expand the neutrino flavor-fluxes in 
terms of the small parameters U e s and ir/4 — 623, and show that there are universal 
first and second order corrections. The second order term can exceed the first order 
term, and so should be included in any analytic study. We also investigate the 
probabilities and ratios after a further expansion around the tribimaximal value of 
sin 2 9i2 = 1/3. In the second part of the paper, we discuss implications of deviations 
of initial flavor ratios from the usually assumed, idealized flavor compositions for pion, 
muon-damped, and neutron beam sources, viz., (y e : i/„ : v T ) = (1:2:0), (0:1: 0), 
and (1:0:0), respectively. We show that even small deviations have significant 
consequences for the observed flavor ratios at Earth. If initial flavor deviations are 
not taken into account in analyses, then false inferences for the values in the PMNS 
matrix elements (angles and phase) may result. 
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1 Introduction 



There has been recently much discussion on the flavor mixing of high energy astrophysical 
neutrinos [1-15]. Neutrino mixing modifies the initial flavor distribution of fluxes $° : 
<3>° : $° in calculable ways. In terrestrial neutrino telescopes such as IceCube [16] or 
KM3Net [17], one can measure neutrino flavor ratios [18] and thereby obtain information 
on the neutrino parameters and/or the sources 1 . There are two essential ingredients to 
such analyses. One is the initial flavor composition which depends on the nature of the 
production process at the source. The other is the neutrino mixing scheme, in particular 
the values of the parameters governing neutrino mixing. The latter plays the role of altering 
the flavor mix from the original due to non-trivial lepton mixing. We will discuss in this 
paper precision issues for both these aspects. 

The current best-fit values as well as the allowed la and 3<r ranges of the oscillation 
parameters are [25]: 

•2/1 n Qo+0- 02 : - 08 

Sin 0i2 = U.,W_ a o2, 0.06 > 

sin 2 2 3 = 0.453:R3, (1) 
\U e3 \ 2 < 0.02 (0.04) . 

The CP phase is completely unknown and has a range between zero and 2n. The current 
information on the mixing parameters [26] therefore suggests that 6*13 and the deviation 
from maximal atmospheric neutrino mixing are small parameters and therefore to ex- 
pand the formulae in terms of them. This is quite useful in order to obtain an analytical 
understanding of basically all phenomenological problems of interest. Furthermore, the 
deviation from the value sin 2 612 = | is also small and thereby a third small expansion 
parameter is introduced. An idealized description for the leptonic mixing, or Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix U, which is nevertheless perfectly compatible with 
all experimental information is tribimaximal mixing [27] 



U ~ C/tbm = 



/ 2 V2 \ 

-1 V2 -Vs 
K -1 V2 Vs ) 



(2) 



The transition probabilities of astrophysical neutrinos are functions of the elements of U, 
whose usual PDG parameterization is given in the Appendix. Expanding these probabilities 
in terms of the small parameters motivated by the observed oscillation phenomenology is 
one of the purposes of this paper. Previously, expansions in terms of \U e s\ and 7r/4 — #23 
have been discussed. In this case, we will show here first that in addition to the known 
universal first order correction [8, 11] there is a second order universal correction which can 
exceed the first order one and has to be included in analytical studies. Furthermore, by 



1 In principle, more exotic neutrino properties such as neutrino decay, Pseudo-Dirac structure, magnetic 
moments, interaction with dark energy, or breakdown of fundamental symmetries could also be studied 
[19-24]. 
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expanding the probabilities in terms of the deviation from sin 2 9\ 2 = | we find extremely 
compact expressions. 

In what concerns the initial flavor mix at the source, there are three simple and - as 
we will make clear in this paper - idealized possibilities: 

• The most conventional one is neutrino emission from purely hadronic processes, such 

as p + p — > 7r T — > yU T + i/^ — > e T + i/J +z/ M + Vp, in which the initial neutrino flavor 
mix is identical to the one in atmospheric neutrinos: 

$0 ; $J ; $0 = 1 ; 2 : . (3) 

If the dominant decay process is p + 7 — > 7r + + X then the flavor mix is still 1:2:0, 
however the initial u e flux is absent. This can in principle be checked by taking 
advantage of the Glashow resonance reaction z/ e + e~ — > W~ [1,3,28], which occurs 
at an incident z/ e energy of 6.3 PeV. 

• The muons may lose energy so that the \>l flux is depleted at the energies of interest. 
This can happen in a variety of ways [14,29-34]; e.g., muons can lose energy in 
strong magnetic fields, or get absorbed in matter. In these cases the initial flavor 
composition is simply 

$0 : $0 : $0 = : 1 : (4) 

without any electron or tau neutrinos. Specific models have been discussed. In 
general, one expects that any pion source will have a transition to a muon-damped 
composition, with the transition energy depending on the source properties. Qualita- 
tively, muon damping occurs when the muon's energy-loss length ~ dE/dx * s shorter 
than its decay length (E/m^) cr , i.e., when dE/dx ?Z cr /m^. 

• The third case is that of sources which emit dominantly neutrons. These neutrons 
originate in the photo-dissociation of heavy nuclei. The decays of the neutrons give 
rise to an initial pure "/5-beam" of I7 e [35], i.e., 

$° : $J : $° = 1 : : , (5) 

with no v e or any muon or tau neutrinos. 

One other expected source of both pion-decay neutrinos and neutron-decay neutrinos, in 
separated energy regions [36], is the GZK nucleonic reaction pcr + 7cmb — ► A + — > n + n + . 
The subsequent pion and muon decays produce neutrinos with energies about 20 times 
(i.e., m n /m N x 1/4) below E GZ k ~ 5 x 10 19 eV, while neutron decay produces 17 e 's with 
energies about a thousand times (i.e., /9-decay Q-value/m n ) below Eqzk- 

We will emphasize in the present work that the above three flavor mixes in Eqs. (3, 
4, 5) are idealizations. Realistically, one should expect deviations from these simple flux 
compositions. These deviations should be taken into account in analyses to avoid incorrect 
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conclusions about the violation/conservation of CP, the octant of #23, or the magnitude of 
\U e s\. We shall outline the need for care with several examples. 

This paper is built up as follows. In Section 2 we discuss an expansion of the mixing 
probabilities and flux ratios up to second order (i) in the small parameters sin 2 #23 — \ and 
\Ue3\1 and eventually, (ii) in the small parameter sin 2 6*12 — §■ In Section 3 we assess the 
validity of the idealized and often used initial flavor compositions given in Eqs. (3, 4, 5). 
We illustrate the effects of impure neutrino mixes with various examples. We sum up and 
conclude in Section 4. 

2 Mixing Probabilities and Flux Ratios up to Second 
Order 

We discuss in this Section approximate formulae for the mixing probabilities and flux 
ratios, and general properties of the relevant flux ratios. For illustration and insight, we 
will expand formulas in terms of small parameters, but we use the exact expressions for all 
plots of flux ratios presented in this paper. 

An expansion up to second order in the small parameters 2 reveals universal correction 
terms [37]. We will expand first in terms of the two parameters sin 2 ^23 — ^ and |C/ e 3|, related 
to breaking of the <-> v T symmetry [38]. To reveal the dependence on these two param- 
eters, it is sometimes useful to fix the solar neutrino mixing with the phenomenologically 
valid relation sin 2 Q\2 — \- 

2.1 Mixing Probabilities 

2.1.1 Expansion in Terms of \U e3 \ and e = ^ — 9 2 s 

The distances to most high energy neutrino sources are quite large compared to oscillation 
lengths \jk = 47r E / Am 2 fe , where Am 2 fc = m 2 — to 2 , is the mass-squared difference and 
E the neutrino energy. Even for energies as high as the GZK-cutoff ~ 5 x 10 19 eV, this 
is so. Consequently, the terms involving the mass-squared differences in the oscillation 
probabilities are effectively averaged out; oscillation probabilities are reduced to mixing 
probabilities. The v a <-> v$ mixing probabilities are 

p a/ 3 = T,\ u ^\ 2 \u^\ 2 ■ (6) 

i 

Starting from an initial flux composition $° : $° : the measurable neutrino flux at 
Earth is given by 

2 The importance of second order terms has also been mentioned in Ref. [12]. 
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If tribimaximal mixing (given in Eq. (2)) is assumed, then the flavor- propagation proba- 
bilities are simply 

10 4 4 



-Ptbm — 73- 

to 



4 7 7 
4 7 7 



(8) 



The flavor ordering for both column and row indices of this symmetric matrix are (e, //, r). 

Neglecting the small parameters n/A — #23 and 13 , the probabilities P a p (and therefore 
the observable flux ratios) are functions solely of the solar neutrino mixing angle. The 
exact expressions for the mixing probabilities (given in the Appendix) are fourth order 
polynomials in the sin^. It is therefore more useful to expand the formulae in terms of 
small parameters and truncate after the quadratic terms. We will first expand in terms of 



U e3 \ and e = 



7T 



23 = 2 -sm 2 e 23 + 0(e 3 ) . 



(9) 



The explicit PDG parametrization of the PMNS matrix is given in the Appendix. Writing 
the result in matrix form yields: 



( P P P 

I 1 ee 1 e(i 1 er 

P P P 

1 e(i 1 (i(l 1 (it 

\ P P P 

\ 1 er 1 (it 1 tt 



(1-26 



12 °12 



Ci o S 



Ci n S 



12 13 12 



1(1 



12 °12 



Ci o S 



12 °12 



Ci 9 S 



Ci 9 S 



12 °12 



C 12 S 12) 2 (1 

1 H _ r 2 2 \ 1 (1 _ 2 „2 

2 c 12 ^12/ 2 1 10 1 



12 °12; 



"12 °12y 



/ 4 -2 -2 \ 



■1(1-2 4,^)11^1 



/ 1 -1 \ 



_ / \ 
2 1 1 +A 1-10 +|A 2 1 -1 
V -2 1 1 J V -1 1 / V -1 1 y 

(10) 

As it must, every row and every column sums to 1. The correction linear in \U e s\ and e is 
given by [8,11] 



A = — sin4# 12 cos 5 |C/ e3 | 
and the correction quadratic in |C/ e3 | and e is 



1 



sin 2 2#i2 e , 



'11) 



■a-2 



sin 2 2^2 cos 2 5 |£/ e3 | 2 + 4 (1 - cos 2 i2 sin 2 12 ) e 2 - sin4#i 2 cos5 |£/ e3 |e (12) 
= (sin 2^i2 cos 5 |C/ e3 | — e cos 26>i 2 ) 2 + 3 e 2 . 

2 

The latter is a new result, as is the observation that the universal second order correction A 
is positive semidefinite. With the current la and (3a) ranges of the oscillation parameters, 

— 2 

one finds the following ranges for A and A : 



-0.104) - 0.043 < A < 0.069 (0.117) , and A 2 < 0.061 (0.179) . 



(13) 



It is seen that the second order correction can exceed the first order correction. Conse- 
quently, the first order correction to the flux ratios alone is not sufficient to accurately 
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describe the phenomenology. The second order correction needs inclusion in analytical 
studies. The reason for the large second order term are the sizable numerical coefficients, 
especially the one in front of e 2 (we have checked that the higher order terms have smaller 
coefficients). 

Explicitly, the individual mixing probabilities are 



p 

1 ee 




(l-2c 2 2 s 2 2 )(l-2|f/ e3 | 2 ) , 






P 




c 2 2 s 2 2 + A + (l-2c 2 2 s 2 2 )|f/ e3 | 2 , 






P 

1 er 




c 2 2 s 2 2 - A + (l-2c 2 2 s 2 2 )|f/ e3 | 2 , 






1 MM 




i(l-44)-A + iA 2 -i(l- 


2 c 12 S 12) 


I^e3| 2 


P 




1 -(l-c 2 l2 s 2 l2 )- 1 -A 2 - 1 -(l-2c 2 l2 


S 12) \ U e3 


|2 

i 


P 

TT 




1(i-c 2 2S 2 2 ) + a + Ia 2 -I(i- 


2 c 12 S 12) 


\Ue3\ 2 



(14) 



This may be compared with the same in the case of exact tribimaximal mixing, given in 

Eq. (8). The explicit terms depending on \U e3 \ 2 are small, because the maximal value of 

(1 — 2c 2 2 s 2 2 ) |?7 e3 | 2 is 0.011 (0.031) at la (3a), which is well below the maximal values of 

|A| and A . In addition, for the initial flavor composition of $° : <3>° : $° = 1 : 2 : 0, 

the terms with explicit \U e3 \ 2 cancel in the fluxes. Consequently, the first and second 

2 

order corrections A and A attain a universal status. P ee and P^ T receive only quadratic 
corrections, where the correction to P ee depends only 3 on # 13 and not on 5 or # 23 . Note 
that the universal second order term shows up only in P^, P^ T and P TT , i.e., in the ji-r 
sector. This will be of interest later when we discuss neutron beam sources, in which these 

2 

probabilities, and therefore A , do not show up in the flux ratios. 

Inserting the tribimaximal value sin 2 9\ 2 = | in the expressions for A and A gives a 
more illustrative result. We define the resulting universal corrections as 

1 

9 

Atbm = - (2cos 2 5|f/ e3 | 2 + 7e 2 - v 7 ^ cos <5 | f/ e3 1 e) . 

2 

If for instance 9i 3 = and |e| is larger than 1/7, then A TBM exceeds |A T bm|- Note that 

2 2 

A = A TB M plus quadratic terms and that A = A TBM plus cubic terms. Hence, the ranges 

2 2 

of A TB M and A TBM are changed little from the ranges of A and A given in Eqs. (13). 

— 2 

It is also interesting to consider the allowed range of the second order correction A in 
the case of a vanishing first order correction A. We find that, at la and (3<r), 

if A = , then A 2 < 0.029 (0.093) . (16) 



^TBM = 7: (V2 cos5\U e3 \ + 4e 



(15) 



3 This remains true for P ee even when oscillations do not average out. 
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In particular if the oscillation parameters lie outside their current la ranges, then the 

mixing probabilities and flux ratios can deviate up to ten percent from their tribimaximal 

values even if the first order correction vanishes. 

Note that in the definition of A in Eq. (11), the factor in front of \ — sin 2 #23 is larger and 

has a smaller range than the factor in front of |C/ e3 | cos 5. To be precise, for the allowed 3a 

range of solar neutrino mixing, | sin A6 12 ranges from 0.12 to 0.21, whereas \ sin 2 29\ 2 ranges 

from 0.38 to 0.48. Consequently, the sensitivity to deviations from maximal atmospheric 

neutrino mixing is better than the sensitivity to deviations from |[/ e3 | = 0. The latter is also 

2 

mitigated by the dependence on the unknown CP phase 5. Comparing A from Eq. (12) 

with the first order parameter A, we note that the factor in front of e 2 in A is larger than 

the one in front of e in A. Also in A , the factor in front of e 2 has a smaller range than the 

one in front of |£/ e3 | 2 : sin 2 29i 2 lies between 0.75 and 0.64, whereas 4(1 — cos 2 6>i 2 sin 2 12 ) 

ranges from 3.06 to 3.25, when the allowed 3a values of 9 12 are inserted. Again, sensitivity 

to deviations from maximal atmospheric neutrino mixing is favored over sensitivity to |£/ e3 |- 

2 

In Fig. 1 we show the minimal and maximal allowed values of A and A as a function of 
the neutrino mixing parameters. There is almost no dependence on 9\ 2 . The strongest 
dependence is on 9 23 . The statistically weak preference for sin 2 # 23 = 0.45 would mean for 
013 = that A ~ 0.02 is positive and A 2 ~ 0.008. 

We show in Fig. 2 the distribution of \U e3 \ cos 5 against sin 2 #23 for several values of A, 
where we allow all mixing angles to vary in their allowed 3a ranges. We also indicate the 
la ranges for sin 2 #23 and for |[/ e3 | cos S if 8 — 0, 7r/4 and 71/ 3. Fig. 3 shows the same for 

2 2 

A . Only if e = |C/ e3 | cos<5 = can A vanish exactly. The first order correction A may 
also vanish, if [8, 11] 

I U e3 1 cos 5 = (sin 2 9 23 - tan 29 12 . (17) 
The tribimaximal value of tan 29i 2 is 2\f2. 

2.1.2 Expansion in terms of |C/ e3 |, e = ^ — 9 2 3, and e' = sin 2 6 12 — | 

We can simplify the expressions even more when we introduce a third small parameter, 
taking advantage of the closeness of sin 2 6*12 to |. Let us define 4 

d = arcsin ^ - 9 12 = JL (I _ sin 2 9 l2 ) + 0(e") . (18) 

The best-fit value of sin 2 12 = 0.32 corresponds to e' = 0.0142. From the previous Sub- 
section it is not difficult to obtain the elements of the flavor-propagation matrix P. They 

4 A related expansion can be found in Rcf. [37]. 
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are [37] 



Pee ^ (10 + 4 A) , P e ,= ^(4-2 A + B) , P eT = L(±-2A-B) , 



^(7 + A-B + C), P^ T = ±(7 + A-C) , P rr = l 



(19) 



where 



|2 



A = 2v / 2e / + 7e /2 -5\U e3 \ 
B = 2 (v^|C/e3| COS5 + 46-4^66' + 7e'|C/e3| costf) , (20) 
C = 4{2\U e3 \ 2 cos 2 5 + 7e 2 - V2e\U e3 \ cos5) . 

2 

Note that C = 9 A TBM , the latter being previously defined in Eq. (15). If 623 is maximal 
and \U e3 \ cos 5 vanishes, then B = C = and only A 7^ 0, unless in addition sin 2 9 12 = |. 

——2 

Just as with A , the correction factor C is positive semi-definite and appears exclusively 
in the /i-t sector. The probability P ee receives contributions only from A, whereas P^ T is 
not corrected by B. Deviations from sin 2 # 12 = | show up mainly in A. 

At la and (3a), we obtain the ranges of these (A, B, C) parameters. The ranges are 

(-0.418) - 0.126 < A < 0.117 (0.284) , 

(-1.956) - 0.778 < B < 1.255 (2.177) , (21) 
< C < 0.465 (1.356). 

The tribimaximal values for these parameters are zero. 

2 

Both expansions, one in terms of (A, A ) and the other in terms of (A, B, C), are 
useful. To disentangle the dependence on 9 12 from the dependences on \U e3 \ and 9 23 — 7r/4, 
the expansion defined by (A, A ) may be more helpful; the (A, A ) expansion has the full 
dependence on 9 12 included in it. On the other hand, the range of e' is smaller than the 
range of e and \U e3 \. Therefore, to a good approximation one may set sin 2 9u to |. Below, 
we concentrate on the dependences on \U e3 \ and e = 7r/4 — #23 • 



2.2 Flavor Ratios 

With the help of the probabilities in Eq. (14) and a given initial flavor composition <3>g : 
$° : $° it is easy to obtain approximate formulae for flux compositions or ratios. We will 
focus here on the ratio of muon neutrinos to the total flux $tot> an d on the ratio of v e to 
v T [18,39]: 

T = i^ and R = ^ . (22) 
$tot $r 

Muon neutrinos with energies > 10 2 GeV can be identified via muons emerging from the 
shower. Electromagnetic showers from u e charged current reactions may be identifiable. At 
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energies > 10 6 GeV, v T can be identified by double-bang or lollipop signatures. Sometimes 
the ratio $ M / ($ e + $ r ) is discussed in the literature. It is trivially related to our ratio T, 
being equal to T/(l — T). 

In the 6.3 PeV energy region, the V e flux becomes easy to measure, due to the enhanced 
rate of the "Glashow resonance" [1,3,28]. The reaction is V e + e~ — > W~ . Thus we also 
define the ratio 



where $e is the V e flux. 

In the rest of this paper, we will focus mostly on the T and R ratios. With the formulas 
we have given in Section 2.1, all other flux ratios are easy to obtain. 



2.2.1 Pion Sources 

Pion-sources present a initial neutrino flux ratios <3>g : $° : $° = 1 : 2 : 0. It holds in 
this case that at Earth, $ e = § $ tot (P ee + 2P efl ), $ M = § $ tot (P eM + 2P Mi ), and $ T = 
| $ tot (P eT + 2P flT ), These simplify to 

pion sources: ($ e : $^ : $ T ) = (l + 2 A : 1 - A + A 2 : 1 - A - A 2 ) . (24) 

Here we have used the (A, A ) expansion. The single (and small) terms containing only 
|£7 e 3| 2 , which are present besides A and A in Eq. (14), drop out of these expressions. If 
one uses the (A, B, C) expansion, which includes e' as well, then one finds that A drops 
out of the ratios, to leave: 

/ B B C B C\ 

pion sources: ($ e : % : $ T ) = [l + - : 1 - - + - : 1 - - - - J . (25) 

These relations illustrate that deviations from the "canonical" 1:1:1 result are of order 
(A, A ), or (^ (B, C)) and can therefore exceed 10%. Another result from these formulae 
is that the ratio 

!^1 + 2A 2 = 1 + ^C (26) 

is always larger than or equal to one. We checked that this is also true for the full expression. 
Hence, there cannot be more v T than at Earth if the initial flux composition is 1 : 2 : 0. 
Reorganizing the results in Eq. (24), we get the ratios of our interest. They are 5 

T~i(l-A + A 2 ) and R ~ 1 + 3 A + A 2 + 3 A 2 . (27) 
o 

Numerically, using the full expressions, T lies between 0.32 and 0.39, while R ranges from 
0.82 to 1.48. Therefore, deviations of more than 15% for T, and almost 50% for R can be 
expected. 

5 The expression for R does not exactly follow from Eq. (14), but is obtained by evaluating and expanding 
the full fraction. The difference is however negligible. 
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The ratio of electron neutrinos to all other neutrino flavors is interesting in that it 
receives no quadratic corrections. The ratio is 

*' I (1 + 2 A). (28) 



$tot 3 

The tribimaximal values of the ratios T, R, and & e /& tot are clearly |, 1, and |, respectively. 

2.2.2 Muon-Damped Sources 

The fluxes at Earth for muon-damped sources are found to be 

TBM 

muon-damped sources: ($ e : $ M : $ T ) = (P efJL : P m : P Mr ) = (4:7:7) 

\2\)) 

— >(4-2A + B:7 + A-B + C:7 + A-C) , 

where we first insert the tribimaximal mixing values, and then show the result with cor- 
rections. From these, we get the ratios of interest: 

T =P»^±(7 + A-B + C)^^-A + ±A 2 , (30) 

P efl 4-2A+B 1/ , 36^2 

where we have inserted on the far right-hand sides, the value sin 2 9 12 — |. The quantities A 

2 2 

and A are therefore the quantities Atbm and A TBM defined earlier in Eq. (15). However, 
as mentioned below Eq. (15), the differences are small. For exact tribimaximal mixing we 
have T = and R = |. Refs. [4] and [13] have proposed to use these sources to probe #23 
and the CP phase, respectively. 

2.2.3 Neutron Beam Sources 

Neutron beam sources have an initial 1:0:0 flavor mix. We find for the flavor decompo- 
sition at Earth, 

TBM 

neutron beam: ($ e : ^ : $ T ) = (P ee : P eM : P eT ) = (5:2:2) 

— > (10 + 4A : 4-2A+B : A-2A-B) . 

Again, we give results for exact tribimaximal mixing, followed by the corrected expression. 
From these, we get the flavor ratios of interest: 

T = P eM ~ ^ ( 4 - 2 A + B) - c ?2 4 + A ~ I + A , 



„ P ee 10 + 4A l-2c 2 s 2 / A \ 5/ 9 A \ ( 32 ) 

P=-^~- 77^ o o 1 + ~- 1 + -A 



18 v y 1Z " 9 

4-2A-P ~ c 2 2 s 2 2 ^ + c 2 2 s\ 2 ) ~ 2 V ' 2 

Tribimaximal values are | and |, respectively. Note that C, or alternatively A 2 , which 
show up only in the /x-t sector, do not appear in neutron beam formulae. The ratio of V e 

TBM ^ 

to the total flux, initially unity, is simply Q = P ee — | — > (10 + 4 A)/18. 
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2.2.4 Summarizing the Flavor Ratios 

We summarize the situation for the flux ratios T and R. Tables 1, 2 and 3 show their 
ranges for the currently allowed 3<x ranges of the oscillation parameters. In the case of 
exact tribimaximal mixing, T is 1/3, 7/18 ~ 0.39, and 2/9 ~ 0.22 for pion, muon-damped 
and neutron sources, respectively. If future neutrino oscillation experiments show that 
deviations from tribimaximal are small, then measurements of T with ~ 10% precision 
would distinguish between pion and muon-damped sources. However, a low value of T ~ 
2/9 would clearly indicate neutron sources 6 . On the other hand, currently allowed nonzero 
values of \U e 3\, e or e' lead to a possible overlap of the ratio T for all sources. This is shown 
in Fig. 4, where we display the allowed ranges of T for the 3<r ranges of the oscillation 
parameters. The range at the very left side of the plot is what is relevant to this discussion 
(the relevance of the rest of the plot concerns flux uncertainties, which we introduce in the 
next Section). 

R, the ratio of electron to tau neutrinos, is for tribimaximal mixing 1, 4/7 and 5/2, 
respectively, for the pion, muon-damped, and neutron sources. Hence, much less precision 
would be required in order to distinguish the different source types. A large value of R 
would indicate a neutron source. However, the correction terms due to nonzero \U e z\, e and 
e' are seen to have rather large pre-factors, so that again overlap can be expected. Results 
with these uncertainties are shown in Fig. 5. 

We show in Fig. 6 the ratio Q of V e to the total flux. Recall that this ratio is important 
for the V e + e~ — > W~ reaction in the ~ 6.3 PeV energy bin. With the Q observable, 
the neutron source ( "cosmic /9-beam" ) is clearly differentiated from the other source types. 
An interesting proposal [28] is that a measurement of Q as the ratio of the resonant 
V e + e~ — > W~ 6.3 PeV bin to any off-resonance bin can differentiate between two nuances 
of the pion source, namely pp and p~y beam dumps. The latter has no I7 e 's at production. 
Fig. 6, in which neutrino luminosities are assumed equal for all source types, reveals that 
pp and p7 are in principle distinguishable. In fact, it has been argued that the neutrino 
luminosity from pp is larger than that from pj by a factor ~ 2.4 (see [28] and references 
therein). This implies a better statistical determination of Q with the pp origin, and an 
additional signal for discriminating pp and pj. On the other hand, the muon-damped 
source also scales with this 2.4 factor, potentially confusing the discrimination between pp 
and p7 sources. Also, it has been noted that the 77 — > reaction may provide some 

17 e 's in a hot pj environment [40]. 

We note that for pion sources the ratios depend very little on #12. For the other sources, 
the dependence is stronger, as (unlike for pion sources) in the limit of e = U e s = the ratios 
depend on 12 . In general, however, for equal deviations from the tribimaximal values, the 
dependence on 6 12 is weaker than the dependences on #23 and \U e s\ cos 5. Nevertheless, we 

6 Alternatively, the same statements apply to the ratio ^^/(^e + $ T ), which is 1/2, 7/11 ~ 0.64 and 
2/7-0.29. 
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note that schematically one may write 

Ratio(pion) = c + C\ A + c 2 A , 
Ratio(muon-damped) = f (9 12 ) + f ± (9 12 ) A + f 2 {9 12 ) A 2 , (33) 
Ratio (neutron) = g (9i 2 ) + ^1(^12) A . 

The zeroth order expression for pion sources does not depend on 9\ 2 . In addition, typically 
we have that the magnitude of cq/c^ 2 is larger than the magnitudes of fo{9i 2 ) / J '1,2(^12) 
and go{9u)/ 91,2(612)- For instance, for the ratio T, we find from the above that |co/ci| = 
\c /c 2 \ = 1, whereas |/ //i| = I/0//2I = 5 (1 - c i2 s i2) - 7 / 18 and \do/9i\ = c? 2 sf 2 ~ 2/9. 
Therefore, the ratios for muon-damped and neutron sources are significantly more sensitive 
to A (and A ), and thus on deviations from vanishing U e3 and maximal 9 23 . This simple 

fact is the reason that flux ratios of muon-damped and neutron beam sources are better 

2 

suited to probe such deviations. Noting further that A and A depend strongly on 9 23 , we 
can expect that neutrino telescopes will be most sensitive to the 9 23 parameter. 

3 Uncertainties in Initial Flavor Composition 

Up to this point we have obtained expressions for the flux ratios in the cases of exact initial 
flavor composition. However, as indicated in the Introduction, one expects on general 
grounds some deviations from the idealized 1:2:0 (pion source), 0:1:0 (muon-damped 
source), or 1 : : (neutron source) flavor ratios. The implications of non-idealized source 
ratios, on the extraction of the neutrino mixing parameters, have not been previously 
studied in detail. In this Section, we will first estimate the amount of "impurity" in the 
initial flux compositions. After that, we will discuss some examples leading to incorrect 
inferences if care is not taken. 

3.1 Estimating Realistic Flavor-Flux Compositions 

When the source of neutrinos is a hadronic beam dump which produces pions and conse- 
quently neutrinos from n — > /i — > e u e , the pion decay chain, naive counting gives a 
flux ratio for v e to of 1 : 2 and no v T . However, the "wrong-helicity" muon polarization 
from pion decay makes the from its decay softer, thus reducing the effective count. 
This effect depends on the injection spectrum. For a canonical spectral index of a — 2, the 
predicted initial flavor ratio is 1 : 1.86 (see the detailed discussion in Ref. [34]). In addition, 
the production and decay of kaons also produces neutrinos. The kaon decay modes which 
produce charged pions in the final state (e.g., Kg — > 2tt, K ± — > 2tt, Kl — > 3n) give rise to 
the same ratio of 1 : 1.86. 

We now extend the analysis of Ref. [34] to include leptonic and semi-leptonic K decays, 

and production and decay of heavy flavors. The appreciable K ± — * /i ± + chain and the 
three-body semi-leptonic decay modes of K ± and especially Kl also produce neutrinos. 
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The K ± — > /i ± + chain has a v e : z/ M ratio of 1 : 2.8 and the Ifx, semi-leptonic decays give 
a ratio of 1 : 0.75. We use a K/ir ratio of 0.15 in our estimates, although the dependence 
on this ratio is rather weak. When all the modes are added to the pion chain neutrinos, 
and all the branching ratios and neutrino energy spectra are taken into account (following 
the methods described in [34]), the final flavor mix ends up being surprisingly close to 
the value 1 : 1.86, which was obtained by including only the pion chain. Finally, at high 
energies, production of the heavy flavors c and b is expected. The semi-leptonic decay of 
these heavy flavors gives rise to so-called "prompt" neutrinos, characterized by the flux 
ratio for v e : of 1 : 1. There is also a small flux of v T from the decay mode D s — > r u T , 
which has a branching ratio of about 6%. Bottom quark decays also produce some u T , but 
6-production is down compared to c-production by a factor of about 30. 

After taking all cross-sections, decaying species, branching ratios, and neutrino energy 
distributions into account, our final estimate of the corrected flavor mix is 

$° : $J : $° = 1.00 : 1.852 : 0.001 . (34) 

Remarkably, the v e : ratio remains very close to the original ratio 1 : 1.86 due to just 
pion decay. Notice that the v T content is a negligible 0.1% at most. The main source of 
uncertainty in our estimate is the spectral index at injection, a. Fig. 3 of Ref. [34] shows 
that the ratio 1 : 1.86 is modified by a few percent to 1 : 1.9 if a = 1.7, and to 1 : 1.8 if 
a = 2.3. 

We quantify the deviations with the single parameter (: 

$°:$J:$° = 1:2(1-C):0, (35) 

and we can expect ( to be ~ 0.1. The ratio 1 : 1.86 for electron to muon neutrinos 
corresponds to ( = 0.07. 

Damped muon sources result when muon energy loss mechanisms are operative. The 
effects of muon energy loss for the flavor mix have been discussed by a number of authors 
[14, 29-34] for a number of sources (such as gamma ray bursts), and models (such as 
Waxman-Bahcall type). From Refs. [14,31,34] one can draw two general conclusions: first, 
the onset of muon-damping happens abruptly in energy. In Ref. [14], the onset is within 
a factor of 2 to 3 of 10 8 GeV. Second, the v e flux may be severely suppressed, but it 
never goes to zero. It never falls below 4% in Ref. [14], and never below 2% in Ref. [34]. 
Consequently, we parameterize the initial muon-damped sources as (77 : 1 : 0), where i] is 
at least a few percent. 

Finally, for the neutron sources a pion pollution of order 10% has been estimated in 
Ref. [4]. So for the muon-damped and neutron sources, we introduce the single parameter 
i] and use for the initial flux ratios 

$° : <3>° : <3>° = i] : 1 : (muon-damped) , and 1 : 77 : (neutron) . (36) 

For typical values of ( and 77, we show in Tables 1, 2 and 3 the ranges of the flux 
ratios under consideration. The oscillation parameters are varied in their currently allowed 
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3cr ranges. Results for pion sources are the most stable. We also give in the Tables the 
values of the flux ratios if the neutrino parameters are fixed to their tribimaximal values. 
One concludes that impure initial flux compositions should be taken into account when 
discussing the prospects of inferring neutrino parameters with flux ratio measurements. 

It is worth mentioning that another possible source of deviations from idealized flux 
ratios is new physics. For example, 3 + 2 sterile neutrino scenarios [41,42] (still allowed 
even after the MiniBooNE results [43]), can cause deviations of the flux ratios of order 10%, 
thereby interfering with the program of inferring deviations from three-flavor tribimaximal 
mixing [42,44]. 

We note that a different parameterization for the initial flavor ratios is given in Ref . [6] . 
There the unit normalization of the sum of the flavor ratios is emphasized by introduc- 
ing two polar angles for the vector on the unit sphere, $ e : $ M : $ r = sin 2 £ cos 2 £ : 
cos 2 £ cos 2 ^ : sin 2 £. Although this parameterization (the angular (^-parameter of Ref. [6] 
is not to be confused with our small (^-parameter defined in Eq. (35)) and ours are equiv- 
alent, ours does have the advantage that the introduced parameters in Eqs. (35) and (36) 
are small, thereby allowing the perturbative expansions we present in the next subsections. 
Also, the focuses in Ref. [6] and in our work are different. The former emphasizes determi- 
nation of the initial flux composition, i.e., determining £ and (, when the neutrino mixing 
parameters are known with sufficient precision. In our work we include the influence of 
uncertainties from both the initial fluxes and the neutrino mixing parameters upon the 
experimental program, with particular attention paid to the consequences for extraction 
of precise neutrino parameters. 

3.2 Parameterized Pion Sources 

Justifiably neglecting the v T contribution, the initial flux composition from pion sources is 
$° : $° : $° = 1 : 2 (1 - () : 0. In the limit of ( = 0, we recover the ratios of Section 2.2. 
The detectable fluxes behave according to 



$ e oc P ee + 2(l-C)P eM =l-2Cc 2 2S 2 2 + 2A(l-C)-2C(l-2c 2 2S 2 2 )|f/ e3 | 2 , 

$ M oc P eM + 2(l-C)P w = l-A-C(l-c 2 2S 2 2 ) + A 2 + 2CA + C(l-2c 2 2S 2 2 )|f/ e3 

$ r oc P er + 2(l-C)P Mr =l-A-C(l-c 2 2S 2 2 )-A 2 + C(l-2c 2 2S 2 2 )|f/ e3 | 2 , 



plus cubic terms in the small parameters. With the alternate (A, B, C) expansion, we 
have 



$ e oc 1 + --2C(4-2A + P) , 

$ M oc l-^ + ( l-2((7 + A-B + C) , 



(37) 



i oc 1 



18 



B 



C 
~9 



2((7 + A-C) . 
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We show in Fig. 4 the allowed range of the ratio T for pion sources as a function of (, while 
Fig. 5 shows the same for R, and Fig. 6 the same for Q. The largest effect is seen in the 
ratio Q, in the case of a p~f neutrino source. 

Impure sources will still lead to deviations from "pure" values of the flux ratios, even 
if neutrino mixing is exactly tribimaximal: 

t 9-7C 1 (. C\ , p 9-4C , ^ C ^ 

Ttbm = 27^l8C " 3 ( and i ^ BM= 9^7C 2il+ 3- (38) 

The deviation is seen to be stronger in -Rtbm- 

Possible nonzero #13 and non-maximal #23 m &y be compensated by the "impurity factor" 
(. The flux ratio $ e / $ to t illustrates this "confusion" . From Eq. (28) we know that for ( = 
the ratio is given by |(1 + 2 A). From Eq. (37), it is (neglecting the very small single terms 
depending on |£4s| 2 ) easily obtained that a nonzero A is compensated by an initial flux 
uncertainty if the following relation holds: 

C = - i 3^ S 2 "~ 9A - (39) 

This value of uncertainty returns the ratio to the value |, even though A ^ 0. The ratio 
of muon and tau neutrino fluxes provides another example. Compensation occurs if 

C = "4 • (40) 
A 

In Fig. 7 we show the distribution of \U e z \ cos 8 against sin 2 623 if the flux ratio T is measured 
to be I and 0.35. We choose a pure source with 1:2:0 and three different impure sources 
motivated by our estimates in Section 3.1. Fig. 8 shows the same for R = 1 and R = 1.1. 
It is obvious that the covered area in parameter space changes considerably when the flux 
composition is varied. Table 1 gives the numerical values of the ranges. Though the upper 
and lower limits are barely modified by varying (, there are easily situations in which 
nonzero ( is dramatic. For example, if T was measured to be | and it was known that 
sin 2 623 = \i then one would infer for an ideal pion source that \Ue^\ cos 5 = 0. If in addition 
\U e s\ was known to be nonzero, then 8 would be 7r/2 and CP-violation would be inferred. 
However, if instead, the true value is ( = 0.1, then \U e3 \ cos 5 ~ 0.06 would hold, and if 
|C/ e3 | was known to be 0.06 then CP may or may not be broken. 

The simple examples given here show that care has to be taken when conclusions about 
neutrino mixing parameters are drawn from flux ratio measurements. 

3.3 Parameterized Muon-Damped Sources 

Here we investigate the impact of the non-idealized initial muon-damped source ratios 
$g : $° : $° = i] : 1 : 0. Figs. 4, 5 and 6 show the minimal and maximal values of the 
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ratios T, R and Q. In general, the dependence on initial flavor deviations is larger here 
than it was with pion sources. The fluxes behave according to 



$ e oc P efl + V P ee = ^(A-2A + B + 2r ] (5 + 2A)) , 

« P w + r / P eM = ^(7 + A- J B + C + r ? (4-2A + J B)) , (41) 

$ r oc P MT + r/P er = ^(7 + ^-^ + ^(4-2^-5)) . 

Table 2 shows the numerical effect of rj ^ for neutrino fluxes from muon-damped sources. 
Comparing with Table 1 for pion sources, one sees that the effect of source uncertainty is 
larger for muon-damped sources compared to pion sources. 
If neutrinos mix tribimaximally, then 

As with pion sources, the parameter dependence is stronger for -Rtbm than Ttbm- 

Fig. 9 displays the distribution of \U e s\ cos S against sin 2 623, taking T to be 7/18 and 
0.42, and for simplicity, fixing sin 2 #12 = |. For the two chosen values of 77 — and 0.1, 
the allowed areas do not meet, which shows again that the sensitivity on impure initial 
fluxes is stronger for muon-damped than for pion sources. Moreover, the dependence on 
the actual value of the ratio T is weaker than the dependence on (. In Fig. 10 we show 
correlated dependences, fixing the ratio R to be 4/7 and 0.7. 

Let us discuss a hypothetical but illustrative example. If all neutrino parameters but 
5 were known exactly, and if T were measured without any uncertainty, then the value 
of cos 5 can be extracted [4, 13]. In Fig. 11 we display this possibility, taking optimistic 
values of the other parameters. We assume a large value of 6 13 and maximal #23 in order to 
maximize the dependence on cos 5. Suppose now that T = 0.397 were measured. Assuming 
that 77 = 0, i.e., a very pure muon-damped source, one would conclude that cos 5 = —0.5, 
thereby inferring leptonic CP- violation. However, if in reality rj = 0.1, then the same 
T M = 0.397 value would instead mean that cos 5 = —1 and CP is conserved. 

Let us present another example using Fig. 12. In Fig. 12 the dependence of T on 
sin 2 023 [5] is shown, for two values of \U e 3\ and for i] = and 0.1. The nonzero rj introduces 
a few percent uncertainty in the flux ratio, especially for sin 2 623 > 0.5. In general, the ratio 
decreases. The curves for different \U e s\ do not meet when 77 = 0; however, for 77 = 0.1 they 
cross each other, thereby destroying the possibility to disentangle the two chosen values of 

l^sl- 
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3.4 Parameterized Neutron Beam Source 



Finally we come to the impact of an initial /3-beam source which may not be pure, i.e. 
$g : <3>° : $° = 1 : rj : 0. In this case the flux ratios at Earth are given by 
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(4-2A + B + V (7 + A-B + C)) , (43) 



The tribimaximal values for oscillation parameters leads to 

14 + 7?? If 3 \ , „ 10 + 4?? 5/ 27 \ 

More generally, we show in Figs. 4, 5 and 6 the minimal and maximal values of the ratios 

T, R and Q (flux composition v e : V e : : v ^ : z/ r : V T = : 1 : : r//2 : : 0) that result 

when the oscillation parameters are allowed to vary over their 3a ranges. As was the case 

with muon-damped sources, deviations from pure flux compositions have more impact on 

neutron sources than on pion sources. And again, the dependence on the actual value of 

the ratio T or R is weaker than the dependence on (. Table 3 confirms these remarks. In 

Figs. 13 and 14 the distributions of \U e3 \ cos 5 against sin 2 6*23 for two characteristic values 

of T (2/9 and 0.26) and R (5/2 and 2) are given. 

2 

Recall that A but not A appears in the flux ratios for neutron sources. This makes 
it possible to give a simple formula for the special case where the effect of nonzero A is 
exactly compensated by a nonzero 77. If 



V = * V ; "^rT' ' l^f , (45) 



A + (1 - 2cf 2 gf 2 ) |j7- e3 | 2 _ a f A , 5 

■ -V 2 s 2 

then the tribimaximal value T = c\ 2 s\ 2 necessarily results. 

We give in Fig. 15 an example of the dependence of T on |C/ e3 | and # 2 3- It can be 
seen that impure initial flux compositions can influence statements on the octant of 6*23- 
For instance, measuring T = 0.2 would rule out sin 2 #23 > 0.5 only if 77 = 0. However, if 
r] = 0.1, then T = 0.2 is compatible with sin 2 #23 = 0.55, and the octant of #23 is different 
from the one inferred if i] = 0. 



4 Summary and Conclusions 

We have considered in this paper neutrino mixing and flux ratios of astrophysical neutrinos. 
We first have expanded the expressions in terms of small parameters e = 7r/4 — #23 and 
\U e 3 1, while leaving Qyi free. The small parameters e and K{f/ e3 } measure the symmetry 
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breaking of <-> v T interchange symmetry. With this expansion, we showed that the first 
and second order corrections which characterize the deviations from /i-r symmetry, A and 
A , appear universally. The universal corrections A and A are given in Eqs. (11, 12). 
Each can take values as large as 0.1. Compact results for the mixing probabilities, in terms 
of A and A , are shown in Eqs. (10, 14). 

The second order term A appears only in the fj,-r sector (therefore it is not relevant for 
flux ratios from neutron beam sources) and is positive semidefinite. Because it can exceed 
the first order term, it is necessary to include it in analytical considerations. It vanishes 
only for e = \U e s\ cos 5 = 0, whereas the first order term can vanish also for nonzero values 
of e and \U e s\. 

In general, if the initial flavor mix is exactly 1:2:0, then neutrino mixing transforms 
these ratios to ($ e : $^ : $ T ) = (1 + 2 A) : (1 - A + A 2 ) : (1 - A - A 2 ). Hence, there are 
always more muon than tau neutrinos upon arrival at Earth. The ratio of muon neutrinos 
to all neutrinos can deviate by more than 15% from the tribimaximal value |, while the 
ratio of electron to tau neutrinos can deviate by up to 50% from the tribimaximal value 1. 

As the solar neutrino mixing parameter sin 2 # 12 is close to the tribimaximal value |, 

we next included e' = arcsin — 9±2 in our set of expansion parameters. With this 
expansion set, three universal corrections A, B and C, defined in [37] and reproduced here 
in Eq. (20), characterize the deviations from tribimaximal mixing. Very concise expressions 
for the mixing probabilities result [37], as seen in Eqs. (19). 

In the second part of this paper, we investigated the purity of initial neutrino-flavor 
ratios expected from three types of astrophysical sources. The initial flavor ratios commonly 
considered in the literature are (i) pion sources (with a complete pion decay chain), having 
initial ratios 1:2:0; (ii) muon-damped sources (initial pions but an incomplete decay 
chain), having initial ratios 0:1:0; and (iii) neutron beam sources (so-called "cosmic f3- 
beams"), with initial neutrino flavor ratios 1:0:0. The idealized flavor ratios of all three 
source types are subject to small but important corrections. We investigated the effects of 
realistic corrections on the extraction of neutrino parameters from measurements of flux 
ratios. We found that the muon-damped and neutron beam sources are more sensitive to 
initial flavor deviations than is the pion source. In addition, the muon-damped and neutron 
beam sources are also more sensitive to deviations of oscillation parameters #23 and C/ e 3 
from 7r/4 and zero, respectively. These sensitivities can be easily seen by considering the 
ratio T of muon neutrinos to the total flux. For sin 2 12 = |, this ratio reads 

|(1-A + A 2 -|C) , pionsource (1 : 2 (1 - C) : 0) , 



*tot ' 18 2 6 



— A + i A — 1 77 , muon-damped source (77 : 1 
I + A + 1 77 , neutron beam source (1 : 77 

where the parameterized initial flavor-ratios are 1 : 2 (1— Q : 0, rj : 1 : and 1 



tively. For pion sources, the zeroth order expression is T = |, and with |A, A 2 , (, rj\ < 0.1, 
deviations can be up to 15% for oscillation-induced effects and of order 1% for impure fla- 
vor mixes. For muon-damped sources, on the other hand, the effect of uncertain oscillation 



0), 
0), 

r] : 0, respec- 
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parameters is up to 30%, and the effect of nonzero i] is more than 5%. The effect on neu- 
tron beam sources may be even more dramatic: the observable T can change by more than 
50% due to deviations from tribimaximality, and by order 10% due to impurities. We have 
also considered the ratio R of electron to tau neutrinos, which in the same sin 2 Qyi = 1/3 
limit reads 



1 + 3 A + A 2 + | , pion source (1 : 2 (1 - () : 0) , 

i? = -| £ ~< f(l + 18A + ^A 2 + yr?), muon-damped source (77 : 1 : 0) , 



(l + § A — H i]j , neutron beam source (1 : 77 : 0) . 
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The magnitude of coefficients reveal that the effects of nonzero A, A , £ or 77 is in general 
stronger on R than on T. To be more quantitative (see Tables 1, 2, 3), oscillation effects 
are up to 50% for pion sources and almost a factor of two for muon-damped and neutron 
sources. Impurities in the initial flux composition of 0.1 lead to deviations in the flux ratios 
of 4%, 20% and 15% for pion, muon-damped and neutron sources, respectively. 

We gave several illustrative examples where the assumption of an idealized, pure initial 
flux ratio may easily (mis) lead to incorrect inferences. Wrong inferences may include the 
octant of 9 2 3, the magnitude of |C/ e 3| and the existence of leptonic CP-violation. We stress 
that in future analyses, the intrinsic flux uncertainty should be taken into account before 
inferences are drawn. 
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A Appendix: Mixing Probabilities 



For the sake of completeness, we give here the explicit forms of the oscillation probabilities. 
The lepton mixing, or PMNS, matrix is parameterized in Particle Data Group format as 



U 



( C12 C13 

-S12 C23 - C12 S23 S13 e l 



S12 C13 

C12 C23 - S12 S23 S13 e* 



V s i2 S23 - C12 c 23 S13 e -C12 S23 - S12 c 23 si 3 e 



S13 e- lS \ 

S23 C13 
C23 C13 J 



(Al) 
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where c^- = cos%, s^- = sin%, and we have omitted the Majorana phases, which are 
irrelevant for neutrino oscillations. Using Eq. (6), one obtains 



e/t 



2 4 { 



C 12 S 12 C 23 



+ (4 + S 



12 y s 13 S 23 



+C12 S12 C 23 S 23 C«5 (C 



12 



'D s 13} , 



where c,5 = cos 5, and 



mm 



1 2 Ci o Ci 



12 ^23 S 23 S 13 



+2 {(4 [(4 + 



4d-l 



(' 



'13 
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'23 



13 °23 \ C 13 S 12 

-2(r 2 r 4 

z \ c 23 c 13 



13 



+ S 



12 



5 23 ; 



c 13 S 23 J C 23 



"12 



r 2 + r 2 s 2 

c 13 ^ c 23 6 12 



.2 
"13 



4 + 1 



+S23 

—2 C12 S12 C23 cs (4 ~~ 4) ( c 
With the help of the identities [11,45] 

Per = Pe^(0 23 -> 23 + 7t/2 Or 9 23 
^ = ^(^23^^23 + ^/2 Or 6 23 



2 \ 2 
S 13 S 23 Sl2 



[4 - 4»]) S is} 



#23 + 3vr/2) , 
#23 + 3tt/2) , 



and the unitary relations 



/IT 



P — P 

J e/i - 1 fi- 



e/j 



/i/i , 



1 - P - P 

x 1 er 1 fir 



P J.9P _|_ P _ 1 — P _P _|_P 



e/t 



/t/t ; 



all other probabilities can be readily obtained. 
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general case 


TBM 


composition 


T = <V$ tot 


R = $ e/$ r 


T = $ M /$ to t 


R = $e/$r 


1:2:0 


0.323^0.389 


0.818-=- 1.476 


0.333 


1.000 


1 : 1.90 : 0.001 


0.321 ^0.386 


0.834-=- 1.493 


0.331 


1.017 


1 : 1.85 : 0.001 


0.321 ^0.384 


0.842-=- 1.503 


0.330 


1.026 


1 : 1.80 : 0.001 


0.320^0.382 


0.850-=- 1.513 


0.329 


1.036 



Table 1: Pion sources: ranges of the ratios T = $ M /$ tot and R = $ e /$ T for the current 
3c ranges of the oscillation parameters. We used different flux compositions, pure 1:2:0 
and different impure cases. The values for tribimaximal mixing are also given. 





general case 


TBM 


composition 


T = V$ tot 


R = $ e /$ T 


T = $ M /$ to t 


R = $e/$r 


0:1:0 


0.33^0.51 


0.34-=- 1.13 


0.39 


0.57 


0.05 : 1 : 


0.33^0.49 


0.40-=- 1.17 


0.38 


0.63 


0.1 : 1 : 


0.33^0.48 


0.46-=- 1.21 


0.37 


0.68 



Table 2: Muon damped sources: ranges of the ratios T = §^/§ tot and R = $ e /$ T for the 
current 3<r ranges of the oscillation parameters. We used different flux compositions, pure 
0:1:0 and different impure cases. The values for tribimaximal mixing are also given. 





general case 


TBM 


composition 


T = V$ tot 


R = $ e /$ r 


T = $ M /$ tot 


R = $ e/$ r 


1:0:0 


0.12^0.35 


1.39 -=-5.35 


0.22 


2.50 


1 : 0.05 : 


0.14^0.35 


1.35 -=-4.73 


0.23 


2.34 


1 : 0.1 : 


0.16^0.35 


1.32 -=-4.26 


0.24 


2.21 



Table 3: Neutron beam sources: ranges of the ratios T = &^/& tot and R = & e /&r for the 
current 3a ranges of the oscillation parameters. We used different flux compositions, pure 
1:0:0 and different impure cases. The values for tribimaximal mixing are also given. 
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Figure 1: The minimal and maximal allowed values of the universal first and second order 

— 2 

parameters A and A as a function of the neutrino mixing parameters. The observables 
not specified in the horizontal axis were varied over their currently allowed la (left) and 
3a (right) ranges. 
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Figure 2: Distribution of |[/ e3 | cos 5 against sin 2 6*23 if A takes certain indicated values. 
Indicated also is the allowed la range of (9 23 and of \U e3 \ cos 5 for (from top to bottom 
above zero) 5 = 0, 5 = it/A and 5 = n/3. The value 5 = ir/2 means \U e %\ cos 5 = 0. 




Figure 3: Distribution of \U e s\ cos 5 against sin 2 #23 if A takes certain indicated values. 
Indicated also is the allowed la range of 6*23 and of \U e s\ cos 5 for (from top to bottom 
above zero) S = 0, 5 = 7r/4 and 5 = tt/3. The value 5 = it/2 means \U e s\ cosS = 0. 
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Figure 4: The minimal and maximal values of the ratio T of muon neutrinos to the total 
flux, obtained by varying oscillation parameters within their 3a ranges. The horizontal 
axis labels deviations (£ or rj) from the idealized flux compositions, parameterized as 1 : 
2 (1 — Q : for pion sources, rj : 1 : for muon-damped sources, and 1 : rj : for neutron 
beam sources. The very left side of the plot is therefore the allowed range for a pure flux 
composition. 
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Figure 5: Same as previous Figure for the ratio R of electron to tau neutrinos. 





4 



0.6 
0.5 
0.4 
» 0.3 
0.2 
0.1 



— py 

— pp 

muon-damped 

neutron 



0.05 



0.1 



0.15 



0.2 



Figure 6: Same as previous Figure for the ratio Q of v e to all neutrinos. 




Figure 7: Distribution of |?7 e3 | cos 5 against sin 2 #23 if the flux ratio 3> M /<&tot is measured 
to be I (left) and 0.35 (right), for different initial flavor compositions. The red circles are 
for pure 1:2:0, the magenta crosses are for 1 : 1.90 : 0.001, the blue squares are for 
1 : 1.85 : 0.001 and the green diamonds are for 1 : 1.80 : 0.001. Indicated also is the allowed 
la range of #23 and of \U e s\ cos 5 for (from top to bottom above zero) 5 = 0, 5 = tt/4 and 
5 = 7r/3. The value 5 = tc/2 means \U e3 \ cos 5 = 0. 




Figure 9: Distribution of \U e3 \ cos 6 against sin 2 2 3 if the flux ratio $ M /$ tot is measured 
to be 7/18 (left) and 0.42 (right); the initial flavor mix is r\ : 1 : 0, and sin 2 # 12 = | is 
assumed. 




Figure 11: Extracted value of cos 5 from an exact measurement of the ratio T of muon neu- 
trinos to the total flux for an initial flux composition r\ : 1 : 0. The oscillation parameters 
are fixed to # 2 3 = 7r /4, sin 2 9 12 = |, and \U e s\ = 0.15. 
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Figure 12: Dependence on sin 2 #23 of the ratio of muon neutrinos to the total flux. The 
red lines are for 0:1:0, while the green lines below are for 0.1 : 1 : 0. We have chosen 
5 = 7r and sin 2 01% — ~. 




Figure 13: Distribution of \U e3 \ cos 5 against sin 2 6*23 if the flux ratio $ M /$ tot is measured 
to be 2/9 (left) and 0.26 (right), for an initial flavor mix of 1 : T] : 0. 




Figure 15: Dependence on # 13 of the ratio of muon neutrinos to the total flux. The red 
lines are for 1:0:0, while the green lines above are for 1 : 0.1 : 0. We have chosen 5 = it 
and sin 2 12 = |. 



